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Abstract. We summarize the results from numerical simulations of mass outflows from AGN. 
We focus on simulations of outflows driven by radiation from large-scale inflows. We discuss 
' the properties of these outflows in the context of the so-called AGN feedback problem. Our 
main conclusion is that this type of outflows are efficient in removing matter but inefficient in 
F ^ 1 removing energy. 
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1. Introduction 



, AGN produce the powerful outflows of the electromagnetic radiation that in tunr can 
drive outflows of matter. These outflows and the central location of AGN in their host 
galaxies imply that AGN can play a key role in determining the physical conditions 
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(e.g., Ciotti & Ostrikerlll997UCiotti et al. 


20091 Sazonov et al.ll2005l ISDringel et al. 2005; 


Fabian et al.l 


2006; iMerloni fc Heinj|2008 


, and references therein). 



CZ3 . 

Outflows from AGN are an important link between the central and outer parts of 
the galaxy because AGN are powered by mass accretion onto a super massive black 
hole (SMBH). Thus the SMBH 'knows' about the galaxy through the accretion flow 
whereas the galaxy knows about the SMBH through the outflows powered by accretion. 
To quantify this connection, let us express the radiation luminosity due to accretion as 

G\ ■ 

OO ■ La = e r c 2 M a , (1.1) 



where we invoke the simplest assumption such that the luminosity is proportional to the 
mass accretion rate (Af a ) and a radiative (or the rest-mass conver sion) efficiency (e r ). 
To estimate M a , one often adopts the analytic formula bv lBondi (1952) who considered 



spherically symmetric accretion from a non-rotating polytropic gas with uniform density 
Poo and sound speed Coo at infinity. Under these assumptions, a steady state solution to 
the equations of mass and momentum conservation exists with a mass accretion rate of 



H ' M b = A47rr^p ooCoo , (1.2) 

where A is a dimensionless p arameter tha t, for the Newtonian potential, depends only 
on the adiabatic index, 7 (cf. iBondi 1952 ). The Bondi radius, ra, is defined as 

r B = -5- 1-3 
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where G is the gravitational constant and M is the mass of the accretor. 

To quantify AGN feedback, one can measure its efficiency in changing the flow of 
mass, momentum, and energy. The mass feedback efficiency e m is defined as the ratio of 
the mass-outflow rate at the outer boundary M out to the mass-inflow rate at the inner 
boundary Mj n , i.e., 

£m = Mout/M n . (1.4) 

Here, we consider only the energy and momentum carried out by matter. Therefore, the 
momentum feedback efficiency (e p ) is defined as the ratio of the total wind momentum 
p w to the total radiation momentum (£ a /c), i.e., 



( P =Pw/ (LJc) 



(1.5) 



Finally, the total energy feedback efficiency et is defined as the ratio between the sum 
of the kinetic power (kinetic energy flux) and thermal energy power (thermal energy 
flux) P t h and the accretion luminosity of the system L a (Eq. [l.lj ). i.e., 



where 



and 



£k = Pk/L a 



£th — Pth/ L a 



(1.6) 



(1.7) 



(1.8) 



It follows that e t = + e t h- 

Several very sophis ticated simulations o f feedback effects were r ecentl y performed, e.g., 



Springel et all (j2005l) (SDH05 hereafter). Il)i Matteo et ail (|2005h . and booth fc Schave 
(2009) (BS09 hereafter). These simulations follow merging galaxies in which many pro- 
cesses were included, for example, star formation, radiative cooling in a complex multi- 
phase medium, BH accretion and feedback. In addition, these local and global processes 
were connected. However, all of this was possible at the cost of crude phcnomcnologi- 
cal realizations of some processes and the spatial resolution at a level larger than tb- 
Consequently, the efficiencies were assumed instead of being computed. 



A main result of these simulations is that the fam ous Mbh_— a relation fe.g-. lFerrarese &: Merritt 
2000; [Gebhardtl Tet al. 20011 iTremainellet al. 2 002) was reproduced. In addition, the BH 
mass is little affected by the details of star formation and supernova feedback. 

This is a great success of the curren t cosmological simulation s models. However, as 
pointed out above (see also others e.g.. iBegelman fc Nathll2005T ) the key feedback pro- 
cesses in the models represent "subgrid" physics. Therefore, in these cosmological and 
galaxy merger simulations AGN feedback cannot be directly related to AGN physics. 

On the other hand, models that aim to provide insights to A GN physics do not include 
galaxy but rather focus on rn or even smaller scales (e.g., Proga 2007t Prqga et al 



20081 IKurosawa" fc Progal lioosl IKurosawa fc Progal l2009at IKurosawa fc Prog al l2009bh . 



Thus they cannot be directly related to AGN feedback on large scales. However, these 
smaller scale simulations can be used directly to measure the feedback efficiencies and in 
turn to quantify the effects that are assume d or parametrized in lar ge scale simulations. 

Here, we summarize the main findings from lKurosawa et al.1 (2009) where we attempted 
to determine if AGN can supply energy in the form and amount required by the cosmo- 
logical and galaxy merger simulations. Our approach was to measure M a , and various 
feedback efficiencies based on direct sim ulations of inflows and out flows in AGN, on 
sub-parsec and parsec scales performed by IKurosawa k, Proga ( 2009bl ) (KP09 hereafter). 
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2. AGN Models in Recent Cosmological Simulations 

Before we present our results, we briefly summarize what is required by the cosmolog- 
ical and galaxy merger simulations. 

2.1. Mass Accretion Rates 



As w e mentioned above, cosmological and galaxy merger simulations (e.g., SDH05; ISiiacki et al 
l2007t BS09), the actual physical process of the mass-accretion onto the BH is not explic- 
itly modeled because of a relatively poor resolution. These simulations rely on a separate 
analytical model to describe the small scale processes. The un resolved accr etion process 
is usually described by a Bondi-Hoyle-Little formulation (e.g., Bondi|[l952t ). 



One can illustrate some important issues related to numerical realization of this process 
by considering a simpler Bondi accretion problem. Then, M a can be written as Eq. (|1.2[) . 
The dimensionless constant A (in Eq. [1.2] ) depends on 7 but is of order of unity. The 
Bondi formula relates M a of a BH located at the center to the gas density and the sound 
speed (or equivalently the temperature) of the gas at a large scale. 

However in cosmological simulations, the Bondi accretion is evaluated as 

Ms — a o (2.1) 

°s 

where p and c s are the density and the sound speed estimated near the BH using the sur- 
rounding smoothed particle hydrodynamics (SPH) gas particles. Note that the expression 
contains "the dimensionless parameter" a which is different from A in Eq. (|1.2p . SDH05 
introduced a parameter to overcome the gap in the scale sizes between the numerical 
resolution and the Bondi accretion regime. In a typical cosmological or galaxy merger 
SPH simulation, the smoothing length (~ 10 3 pc) is much larger than the gravitational 
radius of influence or the Bondi radius, r&, which is ~ 2 pc. If we assume the gas located 
at a large distance, heated by the AGN radiation, is "Comptonized" (T w 2 x 10 7 K), 
the corresponding speed of sound (assuming 7 = 5/3) is relatively high (~ 500 km s -1 ). 

It was found find that a very large factor of a is required for low-mass BHs to grow 
their masses; hence, the problem is not strictly a Bondi accretion problem. Most of the 
AGN feedback model in the cosmological simulations mentioned above assume a constant 
value of a — 100 (see also Table 2 in BS09) except for BS09 who allow a to depend on the 
value of local gas density. The assumption of a very large value of a becomes inadequate 
when the local gas density is higher than that required by the formation of the a cold 
interstellar gas phase, and when the cosmological simulation does resolve the Jean length 
and the Bondi radius (BS09). 

BS09 abandoned the assumption of constant a in Eq. (|2.ip . and introduced the fol- 
lowing parametrization of a. 



1 for iih < nj!j 

(riH/nn)' 3 for n H >n£ 



3 , . . (2-2) 



where «h and are the number density of hydrogen and the critical hydrogen number 
density above which the gas is expected to become multi-phase, and star formation 
is expected t o begin via contraction of gas due to thermo-gravitational instability (cf. 
ISchavell2004h . The critical density is chosen as njj = 0.1cm 3 , i.e., the corresponding 
critical hydrogen density is = 1.7 X 10~ 25 g cm -3 . The best fit models of BS09 to some 
observations (e.g., the Mbh _ ct relation) gives f3 = 2.0. Note that the new parametrization 
of a in Eq. (12.21) provides an additional density dependency of the mass-accretion rate 
in Eq. (|2.ip . In the formulation of BS09, M a steeply depends on the density of the 
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surrounding gas i.e., M cx p 3 , for p > p^ while the formulation of SDH05 and the original 
Bondi accretion model (Eq. |1.2j ) always give a linear dependency i.e., M oc p. In most 
of the AGN accretion models in the cosmological simulations (e.g., SDH05; BS09), the 
highest M a is limited to the Eddington rate, i.e., 

M Edd = 4 " GMBHmp (2.3) 

£ r (7TC 

where m p , e r , ctt and c are the proton mass, the radiative efficiency (the rest mass 
to radiation conversion efficiency), the Thomson cross-section and the speed of light, 
respectively. The Eddington ratio (r) is defined as M a /M^dA- 

To illustrate how M a depends on the de nsity (p) in d ifferent approaches, Fig. 1 (left 
panel) compares results from the models bv lBondil(|l952h . SDH05 and BS09. The mass of 
the BH is assumed as Mbh = 10 s M©. In all three models, the speed of sound c s is set to 
520 km s -2 , which corresponds to that of Comptonized gas temperature T w 2 x 10 7 K 
with 7 = 5/3. In the modified Bondi accretion models of SDH05 and BS09, M a is 
limited by the Eddington rate (Eq. |2.3j ). and the radiative efficiency ef is set to 0.1. The 
dimensionless parameter a — 100 is adopted for the model of SDH05, and /3 = 2.0 is 
adopted in the model of BS09. The figure shows that the mass-accretion rate of SDH05 
is larger than the Bondi accretion rate for p < 10~ 20 gcm~ 3 . On the other hand, M a by 
BS09 is similar to that of the Bondi accretion rate (off by a factor of A = 1/4 in Eq. |1.2j ) 
in the low-density regime ( p < 10 gem ), but it is significantly larger than the 
Bondi accretion rate for the density range of 10 < p < 10~ 20 gcm~ 3 . The densities 
above which the accretion proceeds at the Eddington rate are ~ 9.5 x 10 _23 gcm~ 3 and 
~ 6.3 x 10 _24 gcm -3 for the models of SDH05 and BS09, respectively. However, these 
values change depending on the adopted value of c s . In the Bondi accretion model, M a 
reaches the Eddington rate at much higher density (p ~ 10 -20 gcm~ 3 ). 

2.2. AGN Feedback Models 

Once M a is evaluated, one can estimate the amount of the accretion power that is released 
and deposited. In SDH05 and B S09 (als o in many ot hers, cf. Table 2 in BS09), e r = 0.1 
(e.g., Shakura &: Sunvaevlll973 : see also ISoltan 1982 ) is assumed, and kept constant. A 



higher value of e r (~ 0.2) c an be achieve d in an accretion model with a thin disk and 



a rapidl y rotating BH (e.e.. iThorne 



1974). Recent obs ervations suggest a wide range of 



0. 07 (iMartmez-Sansigre fc Tavlorll2009h . 0.30-0.35 (|Wang et aLll2006a() . On the other 



hand. lCao fc Lil (|2008l ) find e r is relatively low (~ 0.08) for M BH < 10 8 M and relatively 
high (> 0.18) for Mbh ^ 10 9 M©. The exact mechanism coupling the accretion luminosity 
of a BH and the surrounding gas is not well known. Therefore, SDH05 and BS09 simply 
assumed that L a couples only thermally (and isotropically) to the surrounding. Using 
Eq. flUJ , the rate of energy deposition to the surrounding (the AGN feedback rate) in 
SDH05 is 

Et = ef L a = e f e r M B HC 2 (2.4) 
where ei is the efficiency of the AGN energy deposition to the surrounding gas, and is 
a free parameter which is to be constrained by observations. BS09 find the models with 
ef = 0.15 matches observations (e.g., the Magorrian relation and the Mbh — o~ relations) 
very well, and similarly SDH05 find ef = 0.05 matches observations well. 



3. Results from Our Numerical Simulations 

In Kurosawa et al.l ( 20091) . we used our physical two-dimensional (2-D) and time- 
dependent hydrodynamical (HD) simulations of AGN flows to investigate the dependency 
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Figure 1. Left panel: The mass-accretion rates from various models: the Bondi accretion model 
l|Bondilll952l l (solid), the modified Bondi accretion model by SDH05 (dashed) and that by BS09 
(dotted), as described in Eqs. (|l,2|l . (|2,l[l and (|2.2[) . respectively. The mass-accretion rates are 
computed as a function of p at a radius r Q which is much larger than m, i.e., r a ^> j"b. In 
all models, the speed of sound c s at r is set to 520 km s -2 , The mass of the BH is set to 
M B H = 10 8 M Q . While a = 100 is adopted in the model of SDH05, /3 = 2.0 is adopted in the 
model of BS09. In the modified Bondi accretion models of SDH05 and BS09, the mass- accretion 
rates are limited to the Eddington rate (Eq. |2.3| ). Right panel: Comparison of the mass -inflow 
rates found in the HD simulations with those predicted by the Bondi accretion model ((Bondi 
1952) (solid line) and with those adopted by SDH05 and BS09 (dashed line). In the density 
range of the models considered here, M a adopted by SDH05 and BS09 are limited by the 
Eddington rate (Eq. [2.3] ); hence, the line is flat. The mass-inflow rates at the inner boundary 
(circles) and those at the outer boundary (squares) of the computational domain are shown as 
a function of the outer boundary density p . The simulations predict M a that are very similar 
to the Bondi accretion rates, but the models have a less steeper dependency on the density. The 
Bondi mass-accretions rates and those of SDH05 and BS09 are computed for the gas with the 
Comp tonized temperature T = 2 x 10 7 K and with 7 = 5/3. The figure is from iKurosawa et al.l 
|2009h . 

of the BH mass- accretion rate on the surrounding gas density, and to measure the AGN 
feedback efficiencies in converting the accretion luminosity into the outward fluxes of 
energy, momentum and mass. To do this, we simply analyzed the simulations previously 
published in KP09. We note that 2-D simu lations are consistent in m any respects with 
their three-dimensional (3-D) counterparts ( Kurosawa fc Progall2009a ). 



3.1. Mass- Accretion Rates 

The mass- inflow rates at the inner boundary M- ln (n) and those at outer boundary 
M- ln (r a ) from the HD simulations are plotted as a function of the outer boundary density 
p in the right panel of Fig. 1. For a given value of p a , M- ln (r;) and M- ln (r Q ) are not equal 
to each other, but rather Mj n (r\) < M- m (r Q ) because of an outflow. The lowest density 
model is an exception since no outflow is formed in this model. For the higher density 
models, an outflow forms, and not all the material entering from the outer boundary 
reaches the inner boundary. A fraction of gas experiences a strong radiation pressure and 
radiative heating, and the direction of flow changes, forming an outflow. 

The figure also shows M a predicted by the Bondi accretion model (Eq. |1.2j ) and those 
computed from the formulations of SDH05 and BS09 (Eqs. |2.1| and |2.2| ). The outer 
radius is much smaller than that of a typical smoothing scale on a SPH cosmological 
simulation (~ 10 3 pc), and the outer density values used in our simulations are much 
larger than a typical local density at BH in the SPH simulations. The higher density at 
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a 10 pc scale is required to produce an outflow. For example, p Q must be greater than 
2 x 10~ 21 gem -3 , to form an outflow with our system setup. In the density range of the 
models considered here, M a adopted by SDH05 and BS09 are limited by the Eddington 
rate (Eq. [23]); hence, the line is flat (cf. left panel of Fig. 1). 

We do not expect our solutions to reproduce the Bondi density dependency of the 
mass-inflow rate because they include the effects of radiative heating and radiation force. 
However, the right panel of Fig. 1 shows the mass-inflow rates from our models are 
very similar to those of the Bondi rates, i.e., the rates are of the same order of mag- 
nitude. Interestingly our M- ln (r{) and the Bondi mass- accretion rate matches around 
p = 4 x lCP 20 gcm -3 which coincidentally corresponds to V w 1. The accretion rates 
from SDH05 and BS09 are Eddington (the rates corresponding to T = 1) in this density 
range. Therefore their corresponding lines also cross at the same point. 

Fig. 1 (right panel) shows that our models have a weaker dependency of the mass-inflow 
rates on the density than that of the Bondi accretion. The power-law fits of data points 
give the slope q = 0.52 (±0.01) for M- m (n) and q = 0.56 (±0.02) for M in (r Q ), which are 
indeed much smaller than that of the Bondi accretion model, i.e., q = 1 (cf. Eq. |1.2j ). 

3.2. Feedback Efficiencies 

Let us now consider AGN feedback efficiencies in energy, momentum and mass using the 
simulation results, as defined in Eqs. (|1.4|) - (|1.6|) . The models used here show some degree 
of variability (typically ~ 10% level). Therefore, the physical quantities used to compute 
the feedback efficiencies are based on the time averaged values. 

3.2.1. Energy Feedback Efficiency 

Fig. [2] (left panel) shows the energy feedback efficiencies, e t , ek and e t h computed based 
on our models, as a function of the Eddington ratio (r). For systems with relatively low 
Eddington ratio (r< 0.4), e t h > £k- On the other hand, for systems with relatively high 
Eddington ratio (r> 0.6), the kinetic feedback dominates the thermal feedback by a 
factor of ~ 10 to ~ 100. The model with T = 0.2 does not form an outflow, indicating 
an approximate T value below which no outflow forms (with our system setup). 

The energy feedback efficiencies increase as T increases, but the efficiencies saturate 
for r > 1. The total energy feedback efficiency peaks at T « 1 with et ~ 10~ 4 . The 
flattening of the efficiencies for V > 1 is caused by the transition of the inflow-outflow 
morphology to a "disk wind" like configuration for the higher T models (see fig. 4 in 
KP09). Because of the mismatch between the direction in which most of the radiation 
escapes (in polar direction) and the direction in which the most of the accretion occurs in 
the system (the equatorial direction) , the radiatively driven outflows in the disk- wind- like 
configuration cannot increase the outflow efficiency by increasing the accretion luminosity 
or equivalently T. A similar behavior is found in the M ont (^o) - T relation of KP09. 

3.2.2. Mass Feedback Efficiency 

Figure [2] (right panel) shows e m plotted as a function T. The dependency of e m on T is 
very similar to that of the momentum feedback efficiency (not shown here). For T < 1, 
the mass feedback efficiency e m increases as T increases, but it starts to decreases slightly 
beyond T « 1. The efficiency peaks around T = 1 with the maximum efficiency value 
~ 0.4. In other words, about 40% of the total inflowing mass is redirected to an outflow. 
The turn-around of e m values is caused by the transition of the outflow morphology to a 
disk-wind-like configuration (see fig. 4 in KP09) for the larger V models, as in the cases 
for the energy and momentum feedback efficiencies (Sec. 13.2.11) . 

The accretion process in our model is fundamentally different from that of the Bondi 
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Figure 2. Left panel: the efficiencies of converting the BH accretion luminosity L a to the 
rate of energy deposition to the surrounding gas are plotted as a function of the Eddington 
ratio (r). The panel shows the kinetic energy feedback efficiency (circles), the thermal energy 
feedback efficiency e t h (squares) and the total energy feedback efficiency e% = ek + eth (triangles), 
separately (see Eqs. |1.6| . |1.7] . and |1.8| ). The maximum total energy feedback efficiency is 
~ 10 -4 . For the models with relatively low Eddington ratio (T< 0.4), the thermal feedback is 
more efficient than the kinetic feedback (e t h > £k)- For the models with relatively high Eddington 
ratio (r > 0.6), the kinetic feedback is more efficient than the thermal feedback by a factor of 
~ 10 to ~ 100. The model with T — 0.2 does not form an outflow, and the vertical line (dashed) 
at r = 0.2 indicates an approximate F value below which no outflow forms. The flattening of 
the efficiencies at V ~ 1 is caused by the transition of the inflow-outflow morphology to a "disk 
wind" like configuration for the larger T models (see fig. 4 in KP09). Right panel: the mass 
feedback efficiency (e m ) plotted as a function the Eddington ratio (r). The efficiency peaks 
around r = 1 with the maximum efficiency value ~ 0.4, i.e., 40% of the total inflowing mass 
is converted to the outflows. The turn-around of e m values is caused by the transition of the 
outflow morphology to a "disk wind" like configuration for the larger F models (cf. KP09). The 
model with V = 0.2 does not form an outflow, and the vertical line (dashed) at r = 0.2 indicates 
an ap proximate F value below which no outflow forms. The figures are from Kurosawa et al.l 
1120091 ). 

accretion model and those adopted in the cosmological simulations (e.g., SDH05; BS09) 
because in our model an outflow and an inflow can simultaneously be formed while the 
Bondi accretion model can only form an inflow. The model of SDH05 and others can 
form either an inflow or an accretion (but not both simultaneously). 



4. Conclusions 

We performed over thirty number of axisymmctric, timc-dcpcndcnt radiation-hydrodynamical 
simulations of outflows from a slowly rotating (sub-Keplerian) infalling gas driven by the 
energy and pressure of the radiation emitted by the AGN (KP09). These simulations fol- 
low dynamics of gas under the influence of the gravity of the central 10 s M black hole 
on scales from ~ 0.01 to ~ 10 pc. They self-consistently couple the accretion-luminosity 
with the mass inflow rate at the smallest radius (a proxy for M a ). A key feature of the 
simulations is that the radiation field and consequently the gas dynamics are axisym- 
metric, but not spherically symmetric. Therefore, the gas inflow and outflow can occur 
at the same time. 

The relatively large number of simulations allows us to investigate how the results 



depen d on the gas density at the outer radius, p a . In a follow-up paper. Kurosawa et al 



( 20091) . we measure and analyze the energy, momentum, and mass feedback efficiencies 
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of the outflows studied in KP09. We compared our M a -p Q relation with that predicted by 
the Bondi accretion model. For high luminosities comparable to the Eddington limit, the 
power-law fit (M a oc p£) to our models yields q ~ 0.5 instead of q = 1.0 which is predicted 
by the Bondi model. This difference is caused by the outflows which are important 
for the overall mass budget at high luminosities. The maximum momentum and mass 
feedback efficiencies found in our models are ~ 10~ 2 and ~ 10 -1 , respectively. However, 
the outflows are much less important energetically: the thermal and kinetic powers in 
units of the radiative luminosity are ~ 10~ 5 and ~ 10~ 4 , respectively. In addition, the 
efficiencies do not increase monotonically with the accretion luminosity, but rather peak 
around the Eddington limit beyond which a steady state disk-wind-like solution exists. 
Our energy feedback efficiencies are significantly lower than 0.05, which is required in 
some cosmological and galaxy merger simulations. The low feedback efficiencies found in 
our simulations could have significant implications on the mass growth of super massive 
black holes in the early universe. 
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